AD 


TECHNICAL  REPORT  ARCCB-TR-91024 


AN  EXTENSION  OF  MESH 
EQUIDISTRIBUTION  TO  TIME-DEPENDENT 
PARTIAL  DIFFERENTIAL  EQUATIONS 


J.  M.  COYLE 


L 


DVr 

ELF  f 

AUG  2  01991 

a 


JULY  1991 


US  ARMY  ARMAMENT  RESEARCH, 
DEVELOPMENT  AND  ENGINEERING  CENTER 

CLOSE  COMBAT  ARMAMENTS  CENTER 
BENET  LABORATORIES 
WATERVLIET,  N.Y.  12189-4050 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


DISCLAIMER 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated  by  other  authorized 
documents . 

The  use  of  trade  name(s)  and/or  manufacturers )  does  not  constitute 
an  official  indorsement  or  approval. 


DESTRUCTION  NOTICE 

For  classified  documents,  follow  the  procedures  in  DoD  5200. 22-M, 
Industrial  Security  Manual,  Section  11-19  or  DoD  5200. 1-R,  Information 
Security  Program  Regulation,  Chapter  IX. 

For  unclassified,  limited  documents,  destroy  by  any  method  that  wil 
prevent  disclosure  of  contents  or  reconstruction  of  the  document. 

For  unclassified,  unlimited  documents,  destroy  when  the  report  is 
no  longer  needed.  Do  not  return  it  to  the  originator. 


SECURITY  CLASS!  PI  CATION  OF  THIS  PAGE  (Whan  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (and  Subtitle) 

AN  EXTENSION  OF  MESH  EQUIDISTRIBUTION  TO 
TIME-DEPENDENT  PARTIAL  DIFFERENTIAL 

EQUATIONS 

5.  TYPE  OF  REPORT  4  PERIOD  COVERED 

Final 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHOR!*; 

J.M.  Coyle 

8.  CONTRACT  OR  GRANT  NUMBER^*) 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

U.S.  Army  ARDEC 

Benet  Laboratories,  SMCAR-CCB-TL 

Watervliet,  NY  12189-4050 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 

AMCMS  No.  6111. 02. H6 10. 011 
PEON  No.  1A04Z0CANMSC 

11.  CONTROLLING  OFFICE  NAME  AND  AODRESS 

U.S.  Army  ARDEC 

Close  Combat  Armaments  Center 

Picatinny  Arsenal,  NJ  07806-5000 

12.  REPORT  DATE 

July  1991 

13.  NUMBER  OF  PAGES 

12 

<*.  MONITORING  AGENCY  NAME  4  ADDRESS (II  dlltarani  /ram  Controlling  Otllca) 

15.  SECURITY  CLASS,  (of  thl»  rmport) 

UNCLASSIFIED 

\Sm.  DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 

18.  DISTRIBUTION  STATEMENT  (oi  thim  Report)  1 

Approved  for  public  release?  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (ol  tha  abstract  entered  In  Block  20.  II  dlllerent  from  Report) 

ia.  supplementary  notes 

Presented  at  the  Eighth  Army  Conference  on  Applied  Mathematics  and  Computing, 
Cornell  University,  Ithaca,  NY,  19-22  June  1990. 

Published  in  Proceedings  of  the  Conference. 

19.  KEY  WORDS  (Continue  on  rereree  aide  II  necaaeary  mtd  Identity  by  block  number) 

Partial  Differential  Equations 

Finite  Elements 

Mesh  Moving 

Equi distribution 

20.  ABSTRACT  rCoaxtbeja  an  rereree  at*  H  naeaaeery  amt  Identity  by  block  number) 

Mesh  equidistribution  for  static  approximation  problems  is  extended  to  a  domain 
that  admits  time  dependency.  Instead  of  some  quantity  being  equidistributed 
over  a  static  mesh,  the  change  in  this  quantity  is  equidistributed  over  a 
dynamic,  moving  mesh.  Applications  are  for  utilization  in  numerical  methods  to 
solve  time-dependent  partial  differential  equations.  This  mesh  moving  scheme 
is  incorporated  into  a  finite  element  code  which  already  refines  adaptively. 
Comparisons  are  made  for  stationary  and  moving  meshes. 

DO  1473  COITION  or  I  MOV <s  is okoletc  ,  UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWb»n  Data  En I at  ad) 


SECURITY  CL ASSI FIC  ATiON  of  this  PACEflWian  Data  Entmrmd) 


TABLE  OF  CONTENTS 


Page 


INTRODUCTION  .  1 

PROCEDURES  .  3 

RESULTS  .  5 

Case  1  .  6 

Case  2  .  7 

Case  3  .  8 

DISCUSSION  .  9 

CONCLUSIONS  .  9 

REFERENCES  .  11 

LIST  OF  ILLUSTRATIONS 

1.  Mesh  trajectories  for  case  1  .  6 

2.  Mesh  trajectories  for  case  2  .  7 

3.  Mesh  trajectories  for  case  3  .  8 


Accession  For 


NTIS  GRA&I 
DTIC  TAB  □ 

Unannounced  Q 

Just lfi cat  Ion _ 


By - - 

-Distribution/ 


Availability  Codes 
i Avail  and/or 


;oist 


ft' 


Special 


i 


INTRODUCTION 


An  equi distributed  mesh  in  one  space  dimension  is  a  partition  of  a  given 
domain  into  subintervals  such  that  some  given  quantity  is  uniform  over  each 
subinterval.  More  specifically,  given  an  interval  (a,b)  and  a  positive  weight 
function  w(x)  defined  on  (a,b),  then  an  equidistributed  mesh  is  a  partition 

{a  =  x0  <  xi  <  X2  <  . • .  <  xm_i  <  xm  =  b) 

such  that 

*  j  i  p 

f  w(x)dx  =  constant  =  -  f  w(x)dx  ,  j  =  1,2 . M  (1) 

'xj-1  M  'a 

The  usual  application  of  such  a  mesh  is  for  approximating  functional  rela¬ 
tionships  to  a  certain  accuracy  with  a  minimum  number  of  mesh  points  by  choosing 
w(x)  appropriately  (ref  1).  Equidistribution  strategies  have  also  been  used  in 
numerical  methods  for  solving  two-point  boundary  value  problems  (refs  2.3). 

This  is  because  it  has  been  shown  (refs  4,5)  that  the  task  of  selecting  a  mesh 
to  minimize  the  discretization  error  is  asymptotically  equivalent  to  equidis- 
tributing  the  local  discretization  error. 

The  successes  in  the  above  fields  of  functional  approximation  and  numerical 
ordinary  differential  equations  have  led  some  investigators  to  consider  the  use 
of  equidistribution  strategies  for  generating  moving  meshes  in  the  field  of 
numerical  partial  differential  equations  (PDEs)  (refs  6-8).  The  general  frame¬ 
work  is  to  simply  reconsider  Eq.  (1)  with  a  time  dependency.  That  is  to  say. 
the  problem  is  now  to  determine  a  dynamic  mesh 

{a  =  x0  <  x-^t)  <  x2(t)  <  ...  <  xM-i(t)  <  xM  =  b( 
at  time  t  so  that 
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where  the  positive  weight  function  w(x,t)  is  usually  chosen  to  be  a  function  of 
the  solution  of  the  underlying  PDE.  For  example,  w  has  been  chosen  to  be  pro¬ 
portional  to  the  solution's  gradient,  curvature,  and  local  discretization  error. 

When  applying  Eq.  (2)  in  some  numerical  scheme,  most  investigators  move  a 
fixed  number  of  points  to  follow  and  resolve  local  nonuniformities  in  the  solu¬ 
tion.  In  order  to  guarantee  a  certain  accuracy,  they  must  be  sure  that  this 
fixed  number  is  large  enough  to  approximate  the  solution  throughout  the  entire 
spatial  domain  for  the  entire  temporal  "life"  of  the  solution.  Some  see  this  as 
a  limitation  since  the  correct  number  of  fixed  points  necessary  is  not  generally 
known  a  priori. 

Also,  this  moving  mesh  is  not  operating  in  a  vacuum.  It  is  being  used  in 
conjunction  with  some  numerical  solution  procedure.  Since  the  accuracy  of  most 
such  procedures  can  depend  on  the  shape  of  the  soace-time  grid,  sometimes  the 
equidistribution  law  can  be  too  dynamic  and  deform  the  grid  enough  to  introduce 
a  new,  even  larger  source  of  error.  This  can  happen  even  if  the  equidistribu¬ 
tion  law  does  not  demand  much  moving  of  its  own  accord.  If  a  nonequidistributed 
mesh  is  used  as  the  initial  mesh,  then  the  grid  can  deform  drastically  as  the 
moving  mesh  tries  to  relocate  to  the  proper  equidistributed  positions.  In  order 
to  avoid  these  difficulties,  some  investigators  have  abandoned  moving  altogether 
and  developed  local  refinement  methods  (ref  3). 

A  local  refinement  method  is  a  procedure  where  uniform  fine  grids  are  added 
to  coarse  grids  in  regions  where  the  solution  is  not  adequately  resolved. 
Although  they  can  guarantee  a  solution  to  a  prescribed  accuracy,  they  can  be 
costly,  as  they  involve  recomputing  the  solution,  and  they  are  not  as  good  as 
moving  mesh  methods  at  reducing  dispersive  errors  in  the  vicinity  of  wavefronts. 
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The  choice  here  has  been  to  combine  local  refinement  with  mesh  moving  based 
on  equidistribution.  The  purpose  is  twofold.  First,  the  refinement  procedure 
is  incorporated  to  avoid  any  drastic  deformation  of  the  grid  by  the  moving  mesh 
as  well  as  guarantee  a  orescribed  accuracy.  Second,  the  mesh  moving  is  applied 
in  order  to  obtain  as  accurate  a  solution  as  possible  for  any  given  discretiza¬ 
tion  to  delay  the  need  for  refinement  for  as  long  as  possible  and  thus  to  reduce 
the  costs  involved. 

Eouation  (21  as  it  stands,  however,  is  not  easily  partnered  with  a  refine¬ 
ment  scheme.  It  is  too  dependent  on  mesh  position  and  the  number  of  extant  mesh 
points.  Hence,  a  refinement  step  can  disrupt  the  nature  of  the  equidistribution 
and  cause  a  drastic  change  in  the  mesh  dynamics  similar  to  that  caused  by  a 
"bad"  initial  mesh. 

The  attempt  to  overcome  the  difficulty  reported  here  was  to  try  to  extend 
Eq.  (2)  in  such  a  way  that  it  worked  with  the  refinement  procedure  rather  than 
against  it.  It  seemed  that  the  dynamics  of  Eq.  (2)  was  based  on  the  static 
spatial  nature  of  Eq.  (1)  and  an  extension  was  needed  that  incorporated  more  of 
the  time  dependency  of  the  domain  and  solution  process. 

In  the  next  section,  this  extension  of  Eq.  (2)  is  presented  as  well  as  the 
algorithm  for  coupling  the  refinement  and  moving  procedures.  Then,  in  the 
following  section,  results  on  a  series  of  test  cases  are  presented  for  com¬ 
parison.  In  the  Oiscussion  section,  the  characteristics  of  the  extended 
equidistribution  law  are  discussed  in  light  of  the  results  of  the  previous  sec¬ 
tion.  Finally,  in  the  last  section,  some  conclusions  are  presented. 

PROCEDURES 

The  basic  principle  behind  this  extension  of  Eq.  (2)  is  to  try  to 
equidistribute  temporal  properties  as  well  as  spatial.  To  this  end,  consider  a 
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typical  time  interval  of  interest  (0,T)  and  a  discretization 
|0  ■  t0  <  t0  <  ti  <  t2  <  ...  <  tN_j  <  tN  =  T} 

of  that  interval.  Then,  given  any  time  level  tn_i  and  any  mesh 

,  n-1  n-1  n-1  n-1  . 

(a  =  x0  <  x-j  <  X2  <  . . .  <  xM  ®  bf 

at  that  time  level,  requires  the  new  mesh 

(a  =  x0  <  xj  <  X2  <  . . .  <  <  *M  =  bl 

at  the  next  time  level  tn  to  satisfy 


,3  1  ,  rD 

1  n  w(x,tn)dx  =  /  n-1  w(x,  tn_j  )dx  ♦  -  {/a  w(x,tn)dx  -  Ja  wU,^.} 

xj_i  xj-l 


j  =  1.2, 


Note  that  Eq.  (3)  is  not  an  equidistribution  law  in  the  sense  of  Eqs,  (1) 
and  (2).  No  quantity  is  being  held  constant  over  any  subinterval  by  the 
enforcement  of  Eq.  (3).  Rather,  it  is  the  change  in  the  quantity  w  over  each 


new  subinterval 


/  n  n4 
(Xj-1 ) 


that  is  allowed  to  vary  by  a  constant  amount  (which  is  proportional  to  the  total 
change  in  w  from  tn_i  to  tn)  when  compared  to  its  values  over  the  old  subinter- 


.  n-1  n-1. 

( x j_ i , x j  ) 

In  a  sense  then,  it  is  the  time  change  of  this  quantity  that  is  being 
equidistributed. 

Note  also  that  the  relationship  between  old  and  new  meshes  is  not  as 


dramatic  as  in  Eq.  (2).  If 


,  n-ljM 

lxj  }>0 


is  a  "bad”  mesh,  e.g.,  suppose  it  was  readjusted  by  a  refinement  procedure,  then 


Eq.  (3)  simply  requires  that 
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differs  from 


by  a  constant  amount  over  each  subinterval  and  does  not  require  any  drastic 
readjustment  to  a  new  equidistributed  position. 

In  order  to  test  the  performance  of  Eq.  (3)  and  its  postulated  properties, 
it  was  incorporated  into  a  numerical  POE  solver  that  already  implemented  an 
automatic  refinement  strategy.  The  overall  solution  algorithm  is  as  follows: 

1.  Move  mesh  to  next  time  level  according  to  Eq.  (3). 

2.  Solve  POE  using  finite  elements  in  space  and  finite  differences  in 
time. 

3.  Estimate  error  that  occurred  in  the  solution  process. 

4.  If  the  error  is  less  than  or  equal  to  a  prescribed  tolerance  and  the 
time  level  is  less  than  T,  then  go  to  step  1. 

5.  If  the  error  is  greater  than  the  tolerance,  refine  in  either  space  or 
time  or  both,  then  go  to  step  2. 

6.  If  the  time  level  is  greater  than  or  equal  to  T,  then  stop. 

RESULTS 

The  following  POE  was  solved  numerically  for  all  test  cases: 

ut  -  Uxd  +  To  u>  =  256  uxx  '  0  <  *  <  1  -  *  >  0 

u(x,0)  =  tanh  10(x-l)  ,  0  <  x  $  1 

u(0,t)  =  tanh  10(-l+t)  ,  t  ^  0 

u ( 1 ,  t )  =  tanh  lot  ,  t  *  0 

The  exact  solution,  u(x,t)  =  tanh  10(x-l+t),  is  simply  a  wavefront  that  moves 
through  the  spatial  domain  from  right  to  left  as  time  progresses.  Optimally, 


the  mesh  should  try  to  follow  the  front  as  it  moves  across  the  interval  (0,1). 

For  all  cases,  the  L2  norm  of  the  error  was  prescribed  to  be  less  than  a 
tolerance  of  0.01,  an  initial  uniform  mesh  with  11  points  (M  =  11)  and  an  ini¬ 
tial  time  step  of  0.05  (At  =  0.05)  was  input,  and  the  solution  process  was 
allowed  to  proceed  for  75  time  steps. 

Case  1 

For  this  case,  no  movement  was  allowed — only  refinement.  Mesh  trajectories 
are  shown  in  Figure  1.  At  the  end  of  75  time  steps,  N  =  37  and  At  =  0.00218. 


Figure  1.  Mesh  trajectories  for  case  1. 

Horizontal  lines  indicate  a 
temporal  refinement  has  occurred. 
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Mesh  Trajectories  for  Case  1.  Stars  on  x-axis  indicate  the  mesh  input 
before  any  moving  or  refining.  Horizontal  lines  indicate  a  temporal  refinement 
has  occurred  (actual  values  of  At  are  not  shown). 

Case  2 

For  this  case,  movement  was  based  on  the  first  time  derivative  of  the  solu¬ 
tion  (w(x,t)  =  |  ut(x,t)|  ).  Mesh  trajectories  are  shown  in  Figure  2.  At  the  end 
of  75  time  steps,  N  =  40  and  At  =  0.00579. 


Figure  2.  Mesh  trajectories  for  case  2. 

Horizontal  lines  indicate  a 
temporal  refinement  has  occurred. 


Mesh  Trajectories  for  Case  2.  Stars  on  x-axis  indicate  the  mesh  input 
before  any  moving  or  refining.  Horizontal  lines  indicate  a  temporal  refinement 
has  occurred  (actual  values  of  At  are  not  shown). 

Case  3 

For  this  case,  movement  was  based  on  the  second  spatial  derivative  of  the 
solution  (w(x,t)  =  |  uxx(x,t)|  ).  Mesh  trajectories  are  shown  in  Figuie  3.  At 
the  end  of  75  time  steps,  N  =  23  and  At  =  0.00201. 


Figure  3.  Mesh  trajectories  for  case  3. 

Horizontal  lines  indicate  a 
temporal  refinement  has  occurred. 
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Mesh  Trajectories  for  Case  3.  Stars  on  x-axis  indicate  the  mesh  input 


before  any  moving  or  refining.  Horizontal  lines  indicate  a  temporal  refinement 
has  occurred  (actual  values  of  At  are  not  shown). 

DISCUSSION 

Overall,  the  results  are  very  encouraging.  The  mesh  trajectories  flow 
smoothly  with  whatever  solution  characteristic  the  equidistribution  law  is 
based.  This  is  true  even  when  the  initial  mesh  is  unrelated  to  the  equidistri¬ 
bution  rule  and  when  the  refinement  procedure  alters  the  mesh  (see  Figures  2  and 
3).  This  is  exactly  as  desired  and  postulated. 

Furthermore,  it  seems  that  mesh  moving  can  decrease  the  amount  of  refine¬ 
ment  necessary  for  a  given  problem  as  hoped.  This  is  evident  when  comparing 
case  1  with  cases  2  and  3. 

In  case  2,  the  level  of  temporal  refinement  is  less  than  in  case  1  for  the 
same  number  of  time  steps  and  the  same  tolerance  level.  This  is  as  expected 
since  the  temporal  component  of  the  error  is  proportional  to  a  time  derivative 
of  the  solution.  Hence,  movement  based  on  equidistributing  this  error  should 
reduce  the  temporal  refinement  necessary. 

Similarly,  in  case  3,  the  level  of  spatial  refinement  is  less  than  in  case 
1.  Once  again,  this  is  as  expected  since  the  movement  here  is  based  on  a  quan¬ 
tity  proportional  to  the  spatial  component  of  the  error. 

CONCLUSIONS 

Whether  or  not  this  new  moving  scheme  will  develop  into  a  robust  numerical 
procedure  is  still  uncertain.  There  are  still  stability  questions  to  be 
answered  as  well  as  some  implementation  difficulties  not  addressed  here. 
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However,  the  results  presented  here  give  credence  to  the  notion  that  mesh 
moving  and  refinement  schemes  can  be  successfully  combined.  Refinement  proce¬ 
dures  do  not  have  to  interfere  with  mesh  movement,  and  mesh  movement  can  be  per¬ 
formed  to  reduce  the  levels  of  refinement  necessary  to  solve  a  problem  to  a 
given  tolerance. 
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SMCAR-AES ,  BLDG.  321  1 

SMCAR-AET-0 ,  BLDG.  351N  1 

SMCAR-CC  1 

SMCAR-CCP-A  1 

SMCAR-FSA  1 

SMCAR-FSM-E  1 

SMCAR-FSS-D,  BLDG.  94  1 


SMCAR-IMI-I  (STINFO)  BLDG.  59  2 

PICATINNY  ARSENAL.  NJ  07806-5000 

DIRECTOR 

US  ARMY  BALLISTIC  RESEARCH  LABORATORY 
ATTN:  SLCBR-OD-T ,  BLDG.  305  1 

ABERDEEN  PROVING  GROUND,  MD  21005-5066 

DIRECTOR 

US  ARMY  MATERIEL  SYSTEMS  ANALYSIS  ACTV 
ATTN :  AMXSY-MP  1 

ABERDEEN  PROVING  GROUND,  MD  21005-5071 

COMMANOER 
HQ,  AMCCOM 

ATTN:  AMSMC-IMP-L  1 

ROCK  ISLAND,  IL  61299-6000 


NO.  OF 
COPIES 

COMMANDER 

ROCK  ISLAND  ARSENAL 

ATTN:  SMCRI-ENM  1 

ROCK  ISLAND.  IL  61299-5000 

DIRECTOR 

US  ARMY  INDUSTRIAL  BASE  ENGR  ACTV 
ATTN:  AMXIB-P  1 

ROCK  ISLAND.  IL  61299-7260 

COMMANDER 

US  ARMY  TANK-AUTMV  R&D  COMMAND 
ATTN:  AMSTA-DDL  (TECH  LIB)  1 

WARREN,  MI  48397-5000 

COMMANDER 

US  MILITARY  ACADEMY  1 

ATTN:  DEPARTMENT  OF  MECHANICS 
WEST  POINT,  NY  10996-1792 

US  ARMY  MISSILE  COMMAND 
REDSTONE  SCIENTIFIC  INFO  CTR  2 

ATTN:  DOCUMENTS  SECT,  BLDG.  4484 
REDSTONE  ARSENAL,  AL  35898-5241 

COMMANDER 

US  ARMY  FGN  SCIENCE  AND  TECH  CTR 
ATTN:  DRXST-SD  1 

220  7TH  STREET,  N . E . 

CHARLOTTESVILLE,  VA  22901 

COMMANDER 

US  ARMY  LABCOM 

MATERIALS  TECHNOLOGY  LAB 

ATTN:  SLCMT-IML  (TECH  LIB)  2 

WATERTOWN,  MA  02172-0001 


NOTE:  PLEASE  NOTIFY  COMMANOER.  ARMAMENT  RESEARCH,  DEVELOPMENT,  AND  ENGINEERING 
CENTER.  US  ARMY  AMCCOM.  ATTN:  BENET  LABORATORIES,  SMCAR-CCB-TL , 
WATERVLIET,  NY  12189-4050,  OF  ANY  ADDRESS  CHANGES. 


TECHNICAL  REPORT  EXTERNAL  DISTRIBUTION  LIST  (CONT'O) 


NO.  OF  NO.  OF 

COPIES  COPIES 


COMMANOER 

US  ARMY  LABCOM,  ISA 

ATTN:  SLCIS-IM-TL  1 

2800  POWOER  MILL  ROAO 
AOELPHI ,  MO  20783-1145 

COMMANOER 

US  ARMY  RESEARCH  OFFICE 

ATTN:  CHIEF.  IPO  1 

P.O.  BOX  12211 

RESEARCH  TRIANGLE  PARK.  NC  27709-2211 
DIRECTOR 

US  NAVAL  RESEARCH  LAB 
ATTN:  MATERIALS  SCI  &  TECH  DIVISION  1 
CODE  26-27  (DOC  LIB)  1 

WASHINGTON.  D.C.  20375 

DIRECTOR 

US  ARMY  BALLISTIC  RESEARCH  LABORATORY 
ATTN:  SLCBR-IB-M  (DR.  BRUCE  BURNS)  1 

ABERDEEN  PROVING  GROUND,  MO  21005-5066 


COMMANDER 

AIR  FORCE  ARMAMENT  LABORATORY 
ATTN:  AFATL/MN  1 

EGLIN  AFB,  FL  32542-5434 

COMMANOER 

AIR  FORCE  ARMAMENT  LABORATORY 
ATTN:  AFATL/MNF 

EGLIN  AFB.  FL  32542-5434  1 

MIAC/CINDAS 
PUROUE  UNIVERSITY 
2595  YEAGER  ROAD 

WEST  LAFAYETTE.  IN  47905  1 


NOTE:  PLEASE  NOTIFY  COMMANOER,  ARMAMENT  RESEARCH,  DEVELOPMENT,  AND  ENGINEERING 
CENTER.  US  ARMY  AMCCOM,  ATTN:  BENET  LABORATORIES,  SMCAR-CCB-TL . 
WATERVLIET,  NY  12189-4050,  OF  ANY  ADDRESS  CHANGES. 


